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ABSTRACT 

Following the early paper of Goldreich & Julian (1969), polar-cap models have 
usually assumed that the closed sector of a pulsar magnetosphere corotates with the 
neutron star. Recent work by Melrose & Yuen has been a reminder that in an oblique 
rotator, the induction field arising from the time-varying magnetic flux density cannot 
be completely screened. The principal consequence is that the plasma does not corotate 
with the star. Here it is shown that the physics of the polar cap is not changed at the 
altitudes of the radio emission source. But the presence of a plasma drift velocity in 
the corotating frame of reference does provide a mechanism whereby the net charge 
of the star can be maintained within a stable band of values. It also shows directly 
how electron injection and acceleration occur in the outer gap of the magnetosphere. 
It is consistent with radio-loud pulsars in the Fermi LAT catalogue of 7 -emitters all 
having positive polar-cap charge density. 
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1 INTRODUCTION 

The present consensus is that the source of the universal co¬ 
herent radio emission in pulsars is in the open sector of the 
magnetosphere at low altitudes above the polar cap. In this 
region, the B-field is so large that it must be regarded as 
fixed and model-independent, and basic physical considera¬ 
tions can be used to determine the nature of the accelera¬ 
tion field E || and the composition of the plasma. (Through¬ 
out this paper, the parallel and perpendicular directions are 
with respect to the local magnetic flux density B.) The as¬ 
sumption of plasma corotation in the closed sector of the 
magnetosphere has been an essential feature of polar-cap 
models, and follows from the paper of Goldreich & Julian 
(1969) which treated the aligned rotating neutron star with 
no induction component, E 87l “ = 0, in the observer-frame 
electric field. In this case, the electric field generated by the 
rotation is exactly screened by the Goldreich-Julian charge 
density p G j- 

The presence of an induction field in an observable pul¬ 
sar whose magnetic moment is at an angle t[> with the rota¬ 
tion angular velocity has usually not been the subject of 
comment. But two recent papers by Melrose & Yuen (2012, 
2014) have examined the consequences of the induction field 
given directly in the observer frame, in terms of the time- 
dependence of B, by Faraday’s law. They express the field 
as E lnd = E(| nd -1- Ejj 11 * and consider the screening of these 
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components separately. The parallel component can, in prin¬ 
ciple, be screened precisely by a charge density p whilst the 
perpendicular component causes a plasma drift velocity in 
the corotating frame. They note that non-corotation must 
affect the distinction between open and closed sectors of 
the magnetosphere. The consequences of this for polar-cap 
physics and for the outer gap of the magnetosphere are the 
subject of the present brief paper. 

A further problem that has been the subject of little 
comment is the question of how currents flowing from the 
polar cap return so as to maintain the net charge of the star 
within a stable and limited band of values. Polar-cap models 
have ignored this problem, though its treatment is implicit in 
force-free models of the outer magnetosphere. In the compu¬ 
tational aspects of this latter work, the neutron-star radius 
was artificially large, approximately 0.2Rlc, where Rlc is 
the light cylinder radius, and the boundary condition as¬ 
sumed on the neutron-star surface is an equipotential with 
space-charge-limited flow. The force-free solutions specify 
the charge and current densities at all points in the mag¬ 
netosphere. We refer to Bai & Spitkovsky (2010) for some 
informative computed distributions of these quantities. 

The assumption of corotation in the closed magneto¬ 
sphere creates problems in understanding, for example in 
a polar-cap ■ B < 0 pulsar, how the outflow of positive 
charge in the open sector can be compensated by a flux of 
electrons. However, the plasma drift calculated by Melrose 
& Yuen does obviate this difficulty and enables us to see, 
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at least qualitatively, how the necessary charge stability can 
be maintained. 

We refer to Melrose & Yuen (2012, 2014) for a full dis¬ 
cussion of the induction field and related topics but give 
here in Section 2 a very brief summary of the relevant re¬ 
sults. These are considered in relation to the closed magne¬ 
tosphere in Section 3, the polar cap and coherent emission 
processes in Section 4, and the outer gap in section 5. 

Unless otherwise stated, neutron stars with polar-cap 
SI ■ B < 0 are assumed. The reason for this is that the 
polar-cap flux of positive particles has been shown to be 
a non-steady-state plasma, principally of ions and protons. 
The acceleration field E is limited by the reverse flux of 
electrons produced by ion motion through blackbody radi¬ 
ation. The accelerated plasma is relativistic but not ultra- 
relativistic and is an ideal basis for the growth of Langmuir 
modes and turbulence. Pair creation is not significant except 
at very high magnetic flux densities in which case it actually 
suppresses growth of the Langmuir mode. The properties of 
the mode do not otherwise depend on the polar magnetic 
field and so differ little over a wide range of objects, from 
normal pulsars to millisecond pulsars (MSP) and are con¬ 
sistent with the almost universal spectral properties of the 
radio emission. We refer to Jones (2015) and papers cited 
therein for further details of the case for 12 ■ B < 0 pulsars. 
Possible observable consequences of plasma drift in the coro¬ 
tating frame for the SI ■ B > 0 case are considered briefly in 
Section 6. 


2 THE INDUCTION FIELD AND SCREENING 


The induction field E lnd is given directly in the observer 
frame, in terms of the time-dependence of B, by Faraday’s 
law. It is divergence-free, V • E"" i = 0, and so cannot be 
universally screened by the magnetosphere plasma. The par¬ 
allel component can be screened completely, in principle, by 
a time-dependent charge density p = —V • Ej’ ld /(47r) pro¬ 
vided the plasma contains an adequate density of charged 
particles. For the dipole field assumed, this can be expressed 
as, 


P 


~PGJ 0 (cos 6 


COS dm COS Ip) 


9 COS (1 + COS 2 dm) , , 
2(1 + 3COS 2 dm) 2 


(Melrose & Yuen 2014) in which 8 and <j> are observer- 
frame spherical polar coordinates defined with respect to 
the rotation axis and pgjo(v) = £IB/2ttc is the Goldreich- 
Julian charge density defined for 8 = 0 and ip = n, where 
S2 ■ B = QB cos') p. Here rj is the polar coordinate in units of 
the neutron star radius R. Equation (1) is stationary in the 
corotating frame. The magnetic colatitude 8 m is, 


COS 0m = cos ip cos 8 i sin ip sin 8 c.os((p — Lot) , (2) 

in which the appropriate sign is that of the polar-cap SI ■ B. 
The general value of the Goldreich-Julian density is, 

Pgj — t)pgjo (3 cos 8 cos 8m cos if)') . (3) 

The component E|| nd is then completely removed. (The con¬ 
sequent polarization current density given by the continuity 
equation exactly replaces the displacement current compo¬ 
nent <9E|| nd /cdt in Maxwell’s equation.) We shall adopt the 


approximation that p can be treated as a perturbation to 
the aligned corotating state (for small ip, see Melrose & 
Yuen 2014) so that the total screening charge density is 
pa = Pgj + p- 

Melrose & Yuen show that the effect of the compo¬ 
nent E(( ld on the plasma is quite different. In the plane 
locally perpendicular to B, the plasma has a drift veloc¬ 
ity v = cE™ d x B /B 2 relative to the corotating frame. The 
polarization current density in this plane is perpendicular 
to the drift velocity. Owing to the very large values of B 
for a neutron-star plasma, the Alfven velocity is va c 
so that this component of the polarization current density 
is negligible and can give no significant screening of 
The drift velocity given by Melrose & Yuen (2012) for dipole 
field geometry is precisely azimuthal on the magnetic axis 
so that in the vicinity of the polar cap we require only this 
component. Directly from these authors’ equation (10) it is, 


W = ± 


2 QRtj sin ip cos 8 m cos (cp — Sit) 
1 + 3 cos 2 8 m 


( 4 ) 


relative to the corotating frame (Melrose & Yuen 2012) in 
which it is time-independent. Again the correct sign is that 
of SI - B. The degree to which polar-cap physics is affected 
by Vj, is described in Sections 3 and 4. 


3 THE CLOSED MAGNETOSPHERE 

The closed magnetosphere has an atmosphere close to local 
thermodynamic equilibrium. Above that atmosphere, it is 
divided into sectors of positive and negative screening charge 
density separated by a null surface. The neutrality of its 
atmosphere requires that the electron and proton chemical 
potentials satisfy the relation p e (z) — 0) = g P (z) — p p (0) 

at altitudes z within it, and this is achieved through the 
presence of a small non-zero E« which balances the surface 
gravitational acceleration g s , 

e£|| = ~(m p - m e )y, g s < 0. (5) 

The scale height is of the order of 10 _1 cm and the field 
within it is E\\ « 100 V cm -1 . The boundary condition at 
z = 0 allows free transit of electrons, protons or positive ions 
so that electric fields which do not allow an exponentially 
decreasing number density of either electrons or protons at 
the top of the atmosphere are given by eE\Jz) < m e g s or 
ei?||(z) > —m p g a , respectively. 

There are two reasons why such an electric field compo¬ 
nent should exist from time to time to disturb the local ther¬ 
modynamic equilibrium. Firstly, the screening charge den¬ 
sity p a specified by equations (1) and (3) is stationary in the 
corotating frame but is a function of (p and has to be supplied 
by plasma drifting locally perpendicular to B with velocity 
v. This azimuthal variation of p a is shown in Figure 1: the 
extent to which it is significant is an increasing function of 
ip. This takes no account of the additional charge density 
which is required to maintain the open-closed sector bound¬ 
ary condition. The correct charge density on a flux line at 
any instant has to be maintained by transit of electrons or 
protons from the surface atmosphere and a combination of 
drift and motion parallel with B. Also, the flux of protons 
and ions (in an SI ■ B < 0 pulsar) in the open sector leads to 
a rate change of the order of dE^/dt ~ —nv.oPGJoc/2R 2 , in 
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Figure 1. The screening charge density p s is stationary in the 
corotating frame but must be supplied by the drifting plasma. 
It is shown here in units of Pgjo at t = 0 as a function of the 
azimuthal angle o for fixed values of the obliquity 0 and polar 
angle 9. The current flow from the neutron-star surface required 
to maintain the correct charge density is determined by the rate 
of variation of p s as a function of cf>. Each relation is labelled (if), 9) 
in degrees. 

which ito is the open polar-cap radius. This rate of change in 
E« must also be compensated by electron transit at 2 = 0. 
It would be present even if there were no induction field. 

In either case a reverse transfer of protons must also 
be present, but here the very great mass difference becomes 
significant in that the electrons move with relativistic speeds 
at quite low levels of En. Compensation then occurs on any 
flux line at a rate limited only by the velocity of light, as 
described by Melrose & Yuen (2014). In the absence of an 
induction field, there would be a problem in understanding 
how the outward flux of electrons could move perpendicu¬ 
lar to B and so be enabled to cross the light cylinder. The 
presence of the drift velocity v solves this directly. In the 
following Section we attempt to describe the effect of drift 
into the open magnetosphere. 


4 THE POLAR CAP AND OPEN SECTOR OF 
THE MAGNETOSPHERE 

We shall assume that the open magnetosphere divides into 
two sectors as first envisaged by Arons & Scharlemann 
(1979). Its precise definition requires a demarcation between 
those particles that cross the light cylinder and those that do 
not: it depends on knowledge of the E and B fields near the 
light cylinder. Numerical plasma physics, such as the work 
of Bai & Spitkovsky (2010) in the force-free approximation, 
indicates that the cross-sectional area of the open magne¬ 
tosphere is roughly circular at the inner radius of 0.'2riLC 
adopted by them. We refer to flux lines on which the effect 
of the polar-cap acceleration field first described by Mus- 
limov & Tsygan (1992) is not cancelled by flux-line curva¬ 
ture at higher altitudes as the active sector within which the 


energy of outward-accelerated particles is the source of the 
coherent radio emission. For favourable flux-line curvature, 
with a — rf 3 p s (ri) an increasing function of rj, equation ( 6 ) 
shows that further acceleration occurs beyond that of the 
Lense-Thirring effect. The remaining open flux lines will be 
referred to as the return sector. This is obviously not a pre¬ 
cise definition, but those lines on which a zero in fl ■ B exists 
inside the light cylinder are certainly in the return sector. 
This division is shown in Fig. 2 which also shows the drift 
relative to the boundary of the open field lines as it would 
be if there were no inductive field. To be specific, we adopt a 
circular form for the open magnetosphere boundary uo with 
radius uo = 1.6 x 10 4 P ^ 1 ^ 2 cm for a 1.4 Mq star (Harding & 
Muslimov 2001). Above the polar cap and for dipole geome¬ 
try, the open magnetosphere consists of a long narrow tube 
of flux lines whose cross-sectional area given by the open 
sector radius no (??) = 7/ 3 / 2 mo( 1) varies little over lengths of 
the order of the neutron-star radius. Thus if there were no 
division into active and return sectors, the potential at any 
point inside the tube away from the polar-cap surface could 
be expressed as, 

$(u,?7) = n (mq - u 2 ) (P - Ps) • (6) 

This elementary approximation is adequate, particularly as 
the reduction in cross-sectional area caused by the division 
into active and return sectors is uncertain and reduces the 
potential, 4? —¥ /4>, where 1/4 < £ < 1/2 for a semi-circular 
shape. Equation ( 6 ) assumes the condition $ = 0 on the 
open-closed boundary and that it is maintained in the pres¬ 
ence of drift. 

The induction field component ’E‘‘" d causing drift is 
present in the open sector, and its effect on particle motion 
follows at once by considering motion parallel with, and per¬ 
pendicular to, the flux lines. On the basis of equation (4), a 
particle entering an active dipole-field sector at altitude 771 
drifts a distance, 
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( 7 ) 


perpendicular to the flux line it intersected at entry whilst 
moving, inward or outward, with velocity vm to 772- The basic 
unit of length is H — 300P -1 sin-)/) cm for R = 1.2 x 10 6 
cm and dm = c, where P is the rotation period in seconds. 
The velocity in equation (7) is the instantaneous v not its 
average over a rotation period. 

With this evaluation, it is easy to see that drift of pro¬ 
tons or ions into the active open sector from the closed has 
little effect on the polar-cap physics of normal pulsars at 
altitudes less than 77 1 « 10 which contain the source of co¬ 
herent radio emission (Hassall et al 2012). In this case, the 
drift b± at 772 ~ 10 is small compared with the width of the 
active flux tube and the generation of radio frequencies is 
not impaired. This is illustrated by the example shown in 
Fig. 2 for sin 7 / = 1. However, for 772 10, it is possible 

that a fraction of protons accelerated from the neutron-star 
surface will enter the closed sector, but the immediately ob¬ 
servable consequences of this, if any, are unclear. However, 
in the case of millisecond pulsars (MSP), the drift is much 
more significant, although the light-cylinder radius is small, 
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Figure 2. The cross-section of the open magnetosphere at an 
altitude rj = 100 is shown (broken line) divided (equally) between 
the active and return sectors and centred on the magnetic axis 
m as it would be in the absence of drift. Its radius (for dipole 
geometry) is 1.6 X 10 'P” 1 / 2 cm. The displacement of particles 
accelerated from the neutron star surface in the active sector is 
shown by the displaced semi-circle and is 3 X I 0 6 P 1 cm on the 
basis of equation (7). The displacement of particles in the return 
sector which are not accelerated parallel with B is also shown 
and is one half of that in the return sector. The discontinuity 
shown is a result of our approximations and has no consequences 
of significance. 


r/LC ~ 20. Even so, ion and proton acceleration remains 
consistent with radio emission spectra broadly of the uni¬ 
versal form, as is observed (Jenet et al 1998, Kramer et al 
1999, Espinoza et al 2013). The direction of drift is such that 
v$ > 0, that is, super-rotation. It is possible to confirm, di¬ 
rectly from equation (10) of Melrose & Yuen (2012), that its 
direction is identical for both signs of S7 ■ B. 

Drift into the return open sector is a different matter. 
In the absence of an induction field, initial acceleration of a 
test particle arising from the Lense-Thirring effect would be 
nullified in equation ( 6 ) by the unfavourable flux-line cur¬ 
vature, for which < 7 ( 7 ) is a decreasing function of r). Given 
free transit of protons (or ions) at the neutron-star surface, 
the consequence would be stasis unless secondary electron- 
positron plasma exists whose overall charge density could 
adjust so as to accommodate to a null point and an ad¬ 
verse flux-line curvature. But in the absence of pairs, input 
of protons from the positive sector of the closed magneto¬ 
sphere merely contributes to this state and drift continues 
azimuthally through the return open sector into the closed. 

Input of electrons from the closed sector beyond the null 
point on any flux line is followed by outward acceleration 
to the light cylinder. Melrose & Yuen caution that their 
published p and v may be unreliable near the light cylinder, 
but we shall use them to estimate that the total electron 
current passing from the closed sector is of the order of, 

rvLc 

/ dr/Ruo (ri)pGJo 

J Vnull 


= Hcuo(1)pgjo(1) (vie - Vnull) - (8) 

estimated using the instantaneous v. We can see that 
Hr= Rsinipr]^^/ 2 is of the same order as «o(l) so that, 
qualitatively, the electron current given by equation ( 8 ) is 
adequate to compensate for the outward flux of positive 
charge. 


5 THE OUTER GAP 

Although the non-corotation may not be exactly as de¬ 
scribed by Melrose & Yuen, it does have implications for 
the outer-gap model introduced by Cheng, Ho & Ruder- 
man (1986) who first noticed the possible significance of the 
null surface in what we have termed the return sector of the 
magnetospheric charge density distribution. In the SUB < 0 
case, electrons injected beyond it are accelerated outward to 
high energies at the light cylinder and could, in principle, be 
the source of such incoherent radiation as is observed. But 
finding the source of the electrons, or of electron-positron 
pairs, has always been a problem. Single-photon pair pro¬ 
duction by curvature radiation is not usually possible and 
so a further source, photon-photon collisions of curvature 
photons and polar-cap blackbody radiation was later intro¬ 
duced by Cheng, Ruderman & Zhang (2000). 

There has been much work on outer-gap models of 7 - 
emission in the MSP and other pulsars, but particularly in 
the MSP, the above problem remains. The non-corotation 
flux of electrons we have referred to in the final paragraph of 
the previous Section is, in principle, a solution. Electron drift 
into this sector occurs naturally in the SUB < 0 case and the 
integrated flux estimated by equation ( 8 ) is of the same or¬ 
der of magnitude as the Goldreich-Julian active-sector flux 
from the polar cap. It is known that high-energy electrons 
(or positrons) must be present in the vicinity of the light 
cylinder because MSP feature numerously in the Fermi LAT 
catalogue of 7 -emitting pulsars (Abdo et al 2013). But it is 
extremely unlikely that the formation of a dense plasma of 
outward-moving secondary pairs by curvature radiation is 
possible in the outer gap in any model. It must be presumed 
that the outward-moving particles are exclusively electrons 
and that Goldreich-Julian fluxes are sufficient to produce the 
observed 7 -ray spectrum. A further feature of the catalogue 
is the absence of any radio-quiet MSP (defined as a flux 
density S1400 < 30 pjy) which is considered notable even 
though their detection procedures are necessarily different. 
If the outer gap is the source of 7 -emission in the MSP, the 
neutron star must have polar-cap SI • B < 0. This is con¬ 
sistent with the source of the universal radio emission being 
an ion-proton plasma at an S2 ■ B < 0 polar cap (see Jones 
2015). We do not expect radio-frequency emission from an 
SI ■ B > 0 polar cap and in such a neutron star, particle 
drifting into the outer gap region would necessarily be pos¬ 
itively charged, that is, ions or protons, whose acceleration 
would produce negligible incoherent emission. 


6 CONCLUSIONS 

Although the arguments of Melrose & Yuen (2014) are per¬ 
suasive, these authors bear in mind the relative ubiquity of 
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corotation and suggest that it might be enforced by some 
mechanism as yet unknown. But we have discounted this 
possibility and in this paper have taken their results for 
screening charge density and drift velocity as correct, al¬ 
though this cannot be so near the neutron-star surface be¬ 
cause, as they point out, they assume a point-dipole Held or 
near the light cylinder (for example, the denominator term 
in the Goldreich-Julian charge density has been neglected). 
The case in which f2 ■ B < 0 has been chosen for discussion 
because we consider that such neutron stars are the sources 
of the almost universal coherent emission that is observed 
whether in normal pulsars or in MSP (Jones 2013, 2015). 
The fact that screening of the induction field component 
E*" d may be incomplete at radii near the light cylinder is 
likely to have interesting consequences for particle accelera¬ 
tion in that region, as noted Melrose & Yuen (2012). 

In relation to the latter comment, we note that electron 
drift into the region of the return sector beyond the null 
surface at once solves the problem of how the outer gap 
is populated. Observable incoherent emission, particularly 
y-emission is then possible in neutron stars which are not 
capable of supporting pair creation. 

But the polar cap in a magnetosphere that is not in 
strict corotation has, in principle, properties which are un¬ 
changed insofar as radio emission is concerned, certainly if 
the source is at altitudes r/ ~ 10 as favoured by Hassall et al 
(2012). The presence of drift relative to the corotating frame 
also solves, in principle, the return current problem. 

It follows that there appear to be no obvious obser¬ 
vational consequences of non-corotation for radio emission. 
The closed sector is in a dynamic state. Time-dependent 
currents needed to maintain the screening p s and the open- 
closed boundary condition are always present. Particle flow 
inward to the atmosphere is likely to be dissipative but the 
local heating arising from a Goldreich-Julian flux distributed 
over much of the neutron-star surface is unlikely to be ob¬ 
servable. 

The possibility that counter-flowing beams of oppositely 
charged particles might co-exist in a more detailed descrip¬ 
tion of non-corotation is of interest in the case of fl ■ B > 0 
neutron stars. For an ion-proton beam with counter-flowing 
electrons, Langmuir modes exist whose wave-vectors are par¬ 
allel with the ion-proton velocity and so, for this spin sense, 
probably directed outward towards the observer. They re¬ 
quire both beams to have quite small Lorentz factors, but 
the kinematics of the process allows the generation of sub¬ 
stantial radio-frequency power. We refer to Jones (2014) for 
the methods of analysis which lead to this conclusion. The 
observation of unusual or low flux-density sources in the 
Square Kilometre Array (SKA) would be interesting in this 
context. 


ACKNOWLEDGMENTS 

I thank the anonymous referee for some very helpful recom¬ 
mendations which have improved the presentation of this 
work. 


REFERENCES 

Abdo A. A., et al, 2013, ApJ. Suppl., 208:17 

Arons .1., Scharlemann E. T., 1979, ApJ, 231, 854 

Bai X.-N., Spitkovsky A., 2010, ApJ, 715, 1282 

Cheng K. S., Ho C., Ruderman M., 1986, ApJ, 300, 500 

Cheng K. S., Ruderman M., Zhang L., 2000, ApJ, 537, 964 

Espinoza C. M., et al, 2013, MNRAS, 430, 571 

Goldreich P., Julian W. H., 1969, ApJ, 157, 869 

Harding A. K., Muslimov A. G., 2001, ApJ, 556, 987 

Hassall T. E. et al, 2012, A&A, 543, A66 

Jenet F. A., Anderson S. B., Kaspi V. M., Prince T. A., 

Unwin S. C., 1998, ApJ, 498, 365 

Jones P. B., 2013, MNRAS, 435, Lll 

Jones P. B., 2014, MNRAS, 445, 2297 

Jones P. B., 2015, MNRAS, 450, 1420 

Kramer M., Lange C., Lorimer D. R., Backer D. C., Xilouris 

K. M., Jessner A., Wielebinski R., 1999, ApJ, 526, 957 

Melrose D. B., Yuen R., 2012, ApJ, 745-169 

Melrose D. B., Yuen R., 2014, MNRAS, 437, 262 

Muslimov A. G., Tsygan A. I., 1992, MNRAS, 255, 61 

Ruderman M. A., Sutherland P. G., 1975, ApJ, 196, 51 

This paper has been typeset from a TJh'X/ ETj^X file prepared 

by the author. 



